1 



Selection signatures in worldwide Sheep populations - 
Supplementary material 

Maria-Ines Fariello 1 , Bertrand Servin 1 , Gwenola Tosser-Klopp 1 , Rachel Rupp 1 , Carole Moreno 1 , Magali 
SanCristobal 1 , Simon Boitard 3 ' 4 '* 

1 Genetique, Physiologie et Systemes d'Elevage, INRA & ENVT & ENSAT, 
Castanet-Tolosan, France 

3 Genetique Animale et Biologie Integrative, INRA & AgroParisTech, Jouy-en-Josas, 
France 

4 Origine, Structure et Evolution de la Biodiversite, Museum National d'Histoire 
Naturelle & EPHE & CNRS, Paris, France 

* E-mail: Corresponding sboitard@mnhn.fr 



Candidate Gene Populations 



BNC2 


COM 


KITLG 


COM, EMZ 


KIT 


VBS 


EDN3 


MOG 


MC1R 


SUF, SAB, GAL 


ASIP 


MER 


MITF 


VBS, ERS, BOS 



Table SI. List of genes potentially associated with coloration patterns found under selection 
signatures and the most likely candidate populations for selection. 
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Table S2. Parameters for the hapFLK genome scan. K : Number of haplotype clusters used in the 
model for each group, as determined by the fastPHASE cross-validation procedure. [i,a mean and 
standard deviation of the normal distribution used to model the hapFLK neutral distribution 



group DF FLK 

AFR 1 0.97 

ASI 7 6.65 

CEU 3 2.91 

ITA 3 2.93 

NEU 5 4.78 

SWA 3 2.92 

SWE 3 2.96 

Table S3. Theoretical number of degrees of freedom (DF) and observed mean of the FLK statistic 
{FLK) in each group of populations. 
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Figure SI. Projection of animals from the African group on the first 2 principal components. 
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Figure S2. Projection of animals from the North European group on the first 8 principal components. 
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Figure S3. Projection of animals from the Asian group on the first 8 principal components, before 
(top) and after (bottom) correction. One Indian Garolc animal located in the Deccani cluster was 
attributed to this breed. Five Bangladeshi BGE animals clustering away from the rest of the breed were 
removed. Four Changthangi animals, which clustered away from the rest of the breed and appeared 
admixed with the Deccani breed, were removed. One outlier Sumatra animal was also removed. 
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Figure S4. Projection of animals from the Italian group on the first 8 principal components, before 
(top) and after (bottom) correction. Four Altamurana animals located in the Leccese cluster were 
attirubted to the Leccese breed. Similarly, three Lecese animals located in the Altamurana cluster was 
attributed to the Altamurana cluster. 
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Figure S5. Projection of animals from the South West Asian group on the first 8 principal 
components, before (top) and after (bottom) correction. Seven Moghani animals clustering away from 
the rest of the breed, and possibly admixed with the Qezel breed, were removed. 
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Figure S6. Admixture analysis for animals of the Central European group (top), and projection of the 
selected animals on the first 8 principal components (bottom). Three breeds (Swiss Mirror Sheep, Swiss 
Alpine White Sheep and Swiss Black-Brown Mountain Sheep) were reported as admixed breeds in [1] 
and were consequently removed, although they appeared relatively homogeneous in the admixture 
analysis. Two further breeds (East Friesian Brown and Black Headed Mountain) were removed based 
on the admixture analysis. 
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Figure S7. Admixture analysis for animals of the South West European group (top), and projection of 
the selected animals on the first 8 principal components (bottom). For the admixture analysis, a subset 
of 24 animals was sampled at random within each breed in order to obtain balanced sample sizes. Three 
breeds (Castellana, Ojalada and Rasaaragonesa) were removed based on the admixture analysis, as they 
were clearly admixed. Historical records concerning those 3 breeds are ambiguous, some of them 
reporting that they may result from a cross between Merino and Churra animals. The Churra breed 
also looked admixed in this analysis, but this was no longer the case after removing the three other 
problematic breeds. 
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Figure S8. FLK (top) and hapFLK (bottom) p-values in the region of RXFP2. The 
position of RXFP2 is indicated by the vertical gray line. 
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Figure S9. Phylogcnctic tree of the ancestral populations of geographical groups. 
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16 Text SI. FLK and hapFLK genome scans within groups of populations 

17 For each population group, we show, in this order 

is 1. histogram of the observed FLK distribution, and corresponding theoretical x 2 distribution 

19 2. the FLK p-value distribution 

20 3. histogram of the hapFLK distribution and corresponding estimated normal distribution (see details 

21 in Methods) 

22 4. the hapFLK p-value distribution 

23 5. Manhattan plot of the FLK p- values 

24 6. Manhattan plot of the hapFLK p-values 
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FLK genome scan for group AFR 
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hapFLK distribution for group ASI 
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FLK genome scan for group ASI 
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FLK distribution for group CEU 
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hapFLK distribution for group CEU 
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FLK genome scan for group CEU 
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FLK distribution for group ITA 
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hapFLK distribution for group ITA 
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FLK genome scan for group ITA 
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FLK distribution for group NEU 
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hapFLK distribution for group NEU 
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FLK genome scan for group NEU 




FLK distribution for group SWA 
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hapFLK distribution for group SWA 
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P-value distribution for group SWA 
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FLK genome scan for group SWA 
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FLK distribution for group SWE 
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